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STUDY OF TURBOJET COMBUSTOR DYNAMICS USING SWEEP -FREQUENCY DATA 

by John R. Szuch, Francis J. Paulovich, and William M. Bruton 

Lewis Research Center 

SUMMARY 

Tests were conducted to determine the response of the J85 turbojet combustor pres- 
sure to sweep-frequency, sinusoidal oscillations in fuel-spray nozzle pressure. A Fou- 
rier analysis program was used to generate Bode plots of the response characteristics. 

An analog computer simulation of the turbojet combustor was also developed. The 
simulation was used to determine the form of the fuel combustion dynamics required to 
match the experimental data. A dead time second-order lag combination for the fuel 
combustion dynamics was determined. The dead time was relatively insensitive to fuel 
flow, varying between 2. 6 and 3. 2 milliseconds. The second-order lag had a natural 
frequency of 55 hertz and a damping ratio that varied inversely with fuel flow, ranging 
from 0. 7 to 0. 5. The effect of fuel flow on damping ratio was attributed to the primary- 
combustion or burning zone expanding toward the combustor-pressure measuring station. 
The effects of the burning-zone mixture ratio on the response were also determined using 
the simulation. For fixed flows to the combustor, the zero-frequency amplitude ratio for 
the response was a strong function of the burning-zone mixture ratio. 

Transfer functions, obtained from curve fitting the Bode plots of the experimental 
data, were also determined. The response of the combustor pressure to fuel-spray noz- 
zle pressure could be represented by a dead time first-order-lag, second-order-lag 
combination. The dead time was equal to the simulation value for all but the highest fuel 
flow rate case. The first-order lag was related to the mixing of combustion products and 
airflow in the secondary- combustion or mixing zone. The second- order lag had a natural 
frequency of 55 hertz, indicating that the simulated fuel-combustion process dominates 
the closed-loop response that was fit. The damping ratio varied inversely with fuel flow, 
ranging from 0. 9 to 0. 7. 


INTRODUCTION 


The demand for rapid thrust changes in turbojet aircraft requires knowledge of the 
dynamic response of combustor pressure and temperature to disturbances in fuel flow. 
The rapid increase in turbine- inlet temperature must be accomplished while avoiding 
compressor stall, combustor blowout due to excessive fuel-to-air ratios, and excessive 
turbine-blade temperatures. In addition, the analysis of high-frequency interactions be- 
tween engine and inlet dynamics requires an understanding of the combustion dynamics 
involved. 

Evaluations of combustor response to fuel-flow disturbances have often been based 
on transient data. Reference 1 used step changes in fuel flow to study the time response 
of combustor-outlet temperature for comparison with analytical results. Results from 
that study indicated that the response was higher than first order. However, no attempt 
was made to identify the source of the dynamics. References 2 and 3 used similar test 
data to evaluate the effects of various injection methods on performance and blowout 
limits. Reference 4 used step changes in fuel flow to design and evaluate fuel-control 
systems. 

One obvious limitation in trying to evaluate combustor dynamics using transient data 
is the difficulty in measuring small time delays of the order of a millisecond. Refer- 
ence 5 reported results obtained from sinusoidal disturbances in fuel flow. However, no 
phase data were reported, so the observed dynamic response was attributed to lag effects 
rather than to time delays. 

High-response electrohydraulic fuel controls were first developed by NACA in the 
1950's (ref. 6) making the investigation reported in reference 5 possible. Significant im 
provements in these controls have been made in recent years and are reported in refer- 
ences 7 to 9. In addition, improved data analysis techniques have made it possible to 
automatically obtain frequency-response plots when fuel flow is varied sinusoidally with 
frequency being swept linearly or logarithmically. A Fourier analysis program (ref. 10) 
has been developed to reduce the sweep-frequency data yielding amplitude and phase in- 
formation. Accurate phase data at high frequencies (about 100 Hz) allows determination 
of any time delays present in the combustor response. 

Tests were conducted on an afterburning turbojet engine in the Lewis 10- by 10- Foot 
Supersonic Wind Tunnel. Fuel flow was supplied from a high-response electrohydraulic 
fuel valve commanded by a sweep- frequency generator. The Fourier analysis program 
(ref. 10) was implemented on an analog computer. The program yielded Bode plots of 
the response of the combustor pressure to sinusoidal oscillations in the fuel spray-nozzle 
pressure for a range of fuel flow rates. 

For comparison with experimental data an analytical model was developed and imple 
mented on the analog computer. The computer model assumed that the combustor con- 
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sisted of a primary-combustion or burning zone and a secondary-combustion or mixing 
zone. The fuel, together with a portion of the total airflow, was assumed to be completely 
burned in the burning zone. The remaining airflow was assumed to mix with the combus- 
tion products in the mixing zone. The model allowed variations in the burning-zone fuel 
to air ratio. In addition, fuel combustion dynamics could be introduced in the model to 
match the experimental results. 

Transfer functions, obtained from curve-fitting the experimental data, were also 
determined. This information should prove useful in overall system studies where reim- 
ilementation of the entire analog simulation might be undesirable. 


APPARATUS AND PROCEDURE 

A cutaway view of the engine and inlet is shown in figure 1. The engine, a General 
Electric Company model J85-13, is a single-rotor afterburning turbojet engine with an 
ught-stage compressor, an annular combustor, and a two-stage turbine. The engine is 
equipped with a variable-area exhaust nozzle. The inlet is axisymmetric and has mixed- 
•ompression with 60 percent of the supersonic area contraction occurring internally at 
lesign Mach 2. 5. The inlet cowl-lip diameter is 47. 3 centimeters, corresponding to a 
apture flow area of 1760 square centimeters. For this study, only the nonafterburning 
node of operation was used. For all tests, the average free-stream Mach number in 
he wind tunnel was 2. 5 with the inlet- engine system operated at a zero angle of attack, 
'able 1(a) summarizes the test conditions for this study. Fuel flows, ranging from 
. 0517 to 0. 132 kilogram per second, were run with corrected speeds ranging from 83. 1 
3 88 . 3 percent of the design value. 

Figure 2 is a cutaway view of the combustor section. Commercial jet fuel was sup- 
lied to the combustor from a high-response electrohydraulic fuel valve commanded by a 
weep-frequency generator. A description of the valve is found in reference 7. All the 
lel was injected into the main combustor section by means of 12 equally space spray 
ozzles. A cutaway view of one of the nozzles and its flow divider is shown in figure 3. 
he initial flow of fuel passes through the primary tube to the primary-nozzle orifice 
here the spray pattern is formed for combustion. As the pressure across the nozzle 
■creases above 82.7 newtons per square centimeter, a spring-loaded valve in the divider 
msing starts to open. This action allows flow through the secondary tube and orifice to 
ipplement the fuel requirements. Reference 2 discusses the high performance obtained 
ith this type of nozzle. Figure 4 shows a calibration curve for the 12 nozzles obtained 
om cold-flow testing of a single nozzle. 

The pressure measurements of interest for this study were made in the spray-nozzle 
vider housing and in the outer combustor annulus (see fig. 2). Signals from these 
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TABLE I. - STEADY-STATE OPERATING CONDITIONS 


(a) Test measurements 


Test 

Fuel flow, 
kg/sec 

Percent 

corrected speed 

Compressor 
pressure ratio 

Total airflow, 
kg /sec 

Compressor discharge 

Combustor pressure, 
N/cm 3 

Spray-nozzle 


Pressure, 

N/cm 3 

Temperature, 

K 

pressure, 

N/cm 3 

A 

0.052 

86.3 

3.77 

12.6 

30.4 

460 

28.5 

112 

B 

.076 

87. 8 

4. 17 

131 

33.0 

466 

31.8 

123 

C 

.092 

87.4 

4.19 

12.9 

34.0 

479 

32.2 

128 

D 

.106 

83.1 

3.77 

11.6 

30.4 

457 

29.4 

130 

E 

.132 

88.2 

4.63 

12.9 j 

37.0 

485 

35.4 

144 


(b) Simulation values at adjusted fuel-air ratio 


i 

Test 

Burning- zone 
fuel-to-air 
ratio 

Total 

airflow, 

kg/sec 

Burning- zone 
airflow, 
kg/sec 

Mixing- zone 
airflow, 
kg/sec 

Bypass 

airflow, 

kg/sec 

Burning zone 

Mixing zone 


Enthalpy, 

JAg 

Gas 

constant, 

JAg 

Specific 

heat, 

JAg 

Enthalpy, 

JAg 

Gas 

constant, 

JAg 

Specific 

heat, 

JAg 

A 

0.0832 

12.5 i 

0.621 1 

10.6 

1.32 

3.74x10® 

3.60X10 3 

1.63X10 3 

6.53X10 5 

3.45X10 3 

1.04X10 3 

B 

.0918 

12.6 

.830 

10.5 

1.33 

3.41 

3.67 

1.55 

6.97 

3.45 

1.05 

C 

.0932 

12.6 

.993 

10.2 

1.26 

3.38 

3.69 

1.54 

7.55 

3.45 

1.05 

D 

.0863 

11.4 

1.22 

8.53 

1.64 

3.58 

3.62 

1.58 

8.74 

3. 45 

1.07 

E 

.0874 

12.6 

1.51 

9.89 

1.27 

3.55 

3.64 

1.57 

9.16 

3.46 

1.08 
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Figure 3. - Flow divider and spray nozzle. 



Figure 4. - Total fuel flow through 12 spray nozzles as function of spray- 
nozzle pressure drop (based on cold-flow tests of one nozzle). 



transducers were recorded on magnetic tape together with the valve command signal. 

A logarithmic sweep rate of 1 decade per minute was used with frequency ranging 
from 0. 5 to 130 hertz. The frequency response of the fuel valve, when closely coupled 
to the combustor, is flat out to about 100 hertz. For these tests, however, the fuel valve 
was separated from the combustor by about 7 meters of both rigid and flexible lines. 
Resonances in fuel-flow, spray-nozzle pressure, and combustor pressure, which were 
caused by the line dynamics, were observed during the testing. To eliminate these ef- 
fects from the data, the ratio of the response of combustor pressure to disturbances in 
fuel flow to the corresponding response of spray nozzle pressure was used. The strain- 
gage pressure transducers had flat frequency responses (±2 dB) out to about 120 hertz 
for the magnitude of the disturbances used in the test program. 

A Fourier reduction program described in reference 10 was implemented on an ana- 
log computer. The program was used to determine the relation between the combustor- 
pressure oscillations and the spray-nozzle-pressure oscillations. Each previously re- 
corded pressure signal was referenced to the fuel valve command signal to determine, 
as the frequency was varied, the real and imaginary parts of the corresponding transfer 
functions. This information was displayed on an x-y recorder to obtain the amplitude and 
phase angle for each pressure to command signal response. This procedure was re- 
peated for each of the five selected tests. By forming amplitude ratios and subtracting 
the phase angles, the amplitude and phase angle of the combustor-pressure to spray- 
nozzle-pressure ratio transfer function were obtained for each test. A digital reduction 
program has recently been developed to generate the Bode plots directly. Where possi- 
ble, analog-reduced data were compared with digitally reduced data. 


EXPERIMENTAL RESULTS AND ANALYSIS 

The following sections will present the results of the sweep-frequency testing for five 
selected tests. These tests cover a range of fuel flow rates from 0. 0517 to 0. 132 kilo- 
gram per second. The responses of combustor pressure to spray-nozzle pressure, ob- 
tained from the Fourier reduction program, are displayed as Bode plots. 

To match the experimental data and to identify the source of the observed dynamics, 
an analog computer simulation was developed and will be described. Bode plots, obtained 
using the simulation, will be compared with experimental data. 

Transfer functions, obtained from curve fitting of the experimental data will also be 
presented. This information should prove useful in overall system studies where re- 
implementation of the entire analog simulation might be undesirable. 
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TABLE n. - COMPARISON OF EXPERIMENTAL AND SIMULATION DATA FOR J85 COMBUSTOR TEST 


(a) Test A; fuel flow rate, 0. 052 kilogram per second 


[Experimental data normalized at 12 Hz; simulation data normalized at 7 Hz. ] 





Experimental data 



Simulation data 

Frequency, 

Spray-nozzle pres- 

Combustor pres- 

Amplitude 

Phase angle, 

Normalized 

Frequency, 

Amplitude 

Phase angle, 

Normalized 

f, 

sure, P Rn 

sure, P 4 

ratio, 

2 V P sn> 

amplitude 

f, 

ratio, 


amplitude 

Hz 





IV P sJ 

deg 

ratio, 

Hz 

IV^snl 

deg 

ratio, 


Amplitude, 

Phase, 

Amplitude, 

Phase, 



l^sij/ 0 ' 132 




l P 4^snl/ 0132 


N/cm 2 

deg 

N/cm 2 

deg 








7 

2.74 

-50 

0. 374 

-76 

0. 136 

-26 

1.03 

7 

0. 132 

-24 

1.00 

12 

3.05 

-80 

.402 

-119 

. 132 

-39 

1.00 

10 

. 130 

-34 

.985 

17 

3. 01 

“122 

. 370 

-176 

. 123 

-54 

.932 

15 

. 125 

-51 

.947 

22 

2.22 

-154 

.272 

-223 

. 123 

-69 

.932 

20 

. 119 

-66 

.902 

26 

2. 28 

-175 

.258 

-256 

. 113 

-81 

. 856 

25 

. Ill 

-82 

.841 

30 

2.26 

-205 

.238 

-297 

. 105 

-92 

.795 

30 

. 104 

-98 

.788 

34 

2. 33 

-234 

.247 

-334 

. 106 

-100 

. 803 

35 

.0957 

-114 

.725 

38 

2.33 

-265 

.238 

-379 

. 102 

-114 

.773 

40 

.0870 

-128 

.659 

42 

2.05 

-307 

. 192 

-429 

. 0933 

-122 

.707 





46 

1. 37 

-344 

. 112 

-475 

. 0819 

-131 

. 620 

45 

.0782 

-144 

. 592 

51 

.793 

-370 

.0658 

-512 

.0830 

-142 

.629 

50 

.0696 

-158 

.527 

56 i 

. 534 

-386 

.0292 

-533 

. 0547 

-147 

.414 

60 

.0565 

-185 

.428 

74 

.241 

-412 

.0110 

-615 

.0454 

-203 

. 344 

70 

.0464 

-209 

. 352 

78 

.225 

-432 


— 


— 

— 

80 

.0368 

-232 

. 279 

82 

. 192 

-445 


— 


— 

— 

— 


— 


86 

. 165 

-452 


— 


— 

— 

— 



— - 

90 

. 145 

-458 


— 


— 

— 

90 

. 0292 

-253 

. 221 

94 

. 143 

-461 


-740 


-279 

— 

— 


— 

— 

98 

.143 

-466 


-740 


-274 


100 

.0233 

-271 

. 176 


(b) Test B; fuel flow rate, 0. 076 kilogram per second 


Experimental data 


Simulation data 


Frequency, 

Spray-nozzle pres- 

Combustor pres- 

Amplitude 

Phase angle. 

Normalized 

Frequency, 

Amplitude 

Phase angle, 

Nor malized 

f, 

sure, I 


sure, 

P 4 

ratio. 

/P 4 /P sn’ 

amplitude 

f. 

ratio, 

2 V P sn - 

amplitude 

Hz 





IV P snl 

deg 

ratio, 

Hz 

I p 4 /P snl 

deg 

ratio, 


Amplitude, 

2 

Phase, 

Amplitude, 
_ , 2 

Phase, 



l p 4 / p snl / 0 158 




l ^ snl / 0 ' 196 


N cm 

deg 

N/cm 

deg 








7 

6.03 

-34 

0 . 952 

-64 

0.158 

-30 

1.00 

7 

0 . 146 

-25 

1.00 

12 

7.65 

-62 

1. 21 

-106 

. 158 

-44 

1.00 

10 

. 144 

-36 

986 

17 

10.0 

-116 

1. 46 

-172 

. 146 

-56 

. 924 

15 

. 137 

-53 

. 938 

22 

5.48 

-148 

. 752 

-226 

. 137 

-78 

. 867 

20 

. 130 

-70 

. 890 

26 

5 . 34 

-176 

.731 

-269 

. 137 

-93 

. 867 

25 

. 120 

-86 

. 822 

30 

4.77 

-199 

. 614 

-306 

. 129 

-107 

. 816 

30 

. 113 

-102 

. 774 

34 

4.52 

-220 

. 545 

-341 

. 120 

-121 

.759 

35 

. 104 

-118 

.712 

38 

4.52 

-248 

.530 

-380 

. 117 

-132 

.740 

40 

.0956 

-134 

.655 

42 

4. 11 

-282 

. 457 

-426 

. Ill 

-144 

.702 

--- 


— 

— 

46 

3.01 

-316 

. 300 

-472 

. 0995 

-156 

.630 

45 

.0871 

-149 

.596 

51 

1.97 

-340 

. 165 

-512 

.0836 

-172 

. 529 

50 

.0776 

-164 

.531 

56 

1.37 

-354 

. 104 

-541 

.0759 

-187 

.480 

60 

.0631 

-192 

.432 

69 

.658 

-368 

.0421 

-587 

.0639 

-219 

.404 

70 

.0502 

-218 

. 344 

74 

. 814 

-381 

.0457 

-606 

.0562 

-225 

. 356 

--- 


— 

— 

78 

.724 

-397 

.0432 

-621 

.0596 

-224 

. 377 

80 

.0393 

-242 

.269 

82 

.647 

-409 

.0208 

-647 

.0322 

-238 

.204 

— 


— 

— 

86 

.603 

-416 

.0175 

-692 

.0290 

-276 

. 184 

--- 


— 

— 

90 

.603 

-422 

.0183 

-690 

.0303 

-268 

. 192 

90 

. 0309 

-263 

.212 

— 

— 

— 


— 


— 

— 

100 

.0240 

-283 

. 164 
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TABLE n. - Continued. COMPARISON OF EXPERIMENTAL AND SIMULATION DATA FOR J85 COMBUSTOR TEST 
(c) Test C; fuel flow rate, 0. 092 kilogram per second 
[ Experimental data normalized at 12 Hz; simulation data normalized at 7 Hz. ] 

Simulation data 


Experimental data 


Frequency, 

Spray-nozzle pres- 

Combustor pres- 

Amplitude 

Phase angle, 

Normalized 

Frequency, 

Amplitude 

Phase angle, 

Normalized 

f, 

Hz 

sure. 

P sn 

sure, 

P 4 

ratio, 

^V P sn> 

amplitude 

f, 

ratio, 

ivvi 

''V P sn- 

amplitude 

Amplitude, 

N/cm 2 

Phase, 

deg 

Amplitude, 

N/cm 2 

Phase, 

deg 


deg 

ratio, 

l^^snl/ 0 - 126 

Hz 

deg 

ratio, 

Ify^sJ/ 0 - 126 

7 

4.38 

-50 

0. 548 

-80 

0. 125 

-30 

0. 992 

7 

0. 126 

-25 

1.00 

12 

4. 80 

-82 

.603 

-127 

. 126 

-45 

1.00 

10 

. 124 

-35 

.884 

17 

4. 94 

-119 

.567 

-182 

. 115 

-63 

.913 

15 

. 119 

-52 

. 944 

22 

4. 11 

-147 

.438 

-230 

. 107 

-83 

. 849 

20 

. 114 

-69 

.905 

26 

4. 11 

-177 

.430 

-275 

. 104 

-98 

. 825 

25 

. 106 

-86 

. 841 

30 

3.92 

-204 

. 384 

-315 

.0979 

-111 

.777 

30 

. 100 

-102 

.794 

34 

3.70 

-232 

. 353 

-354 

.0953 

-122 

.756 

35 

.0930 

-118 

.738 

38 

3.54 

-261 

.311 

-394 

.0879 

-133 

.698 

40 

.0856 

-134 

.679 

42 

3.01 

-292 

.249 

-434 

.0826 

-142 

.656 

-- 


— 

— 

46 

2.19 

-322 

. 146 

-475 

.0667 

-153 

. 529 

45 

.0779 

-149 

. 618 

51 

1.51 

-345 

.0731 

-513 

.0484 

-168 

. 384 

50 

.0697 

-165 

.553 

78 

.559 

-416 


-640 


-224 


60 

70 

.0563 

.0445 

-193 

-220 

.447 

.353 




82 

.521 

-428 

-680 

-252 


.0344 

-244 

.273 





86 

90 

.510 

.521 

-437 

-447 

-694 

-722 

-257 

-275 







90 

.0267 

-266 

.212 




94 

.548 

-456 

-718 

-262 









.0206 

-286 

.163 

” 









(d) Test D; fuel flow rate, 0. 106 kilogram per second 







Experimental data 




Simulation data 


Frequency, 

Spray- nozzle pres- 

Combustor pres- 

Amplitude 

Phase angle, 

Normalized 

Frequency, 

Amplitude 

Phase angle. 

Normalized 

f, 

Hz 

sure, 

p s „ 

sure, 

P 4 

ratio, 

IV P snl 

deg 

amplitude 

ratio, 

f, 

Hz 

ratio, 

IV P snl 

-V^n' 

deg 

amplitude 

ratio. 






Amplitude, 

N/cm 2 

Phase, 

deg 

Amplitude, 

N/cm 2 

Phase, 

deg 


IVpJ/ 0 - 146 



iP^snl/ 0 - 172 

7 

4.52 

-36 

0.640 

-65 

0.141 

-29 

0.966 

7 

0.172 

-24 

1.00 

12 

5.51 

-64 

. 807 

-104 

. 146 

-40 

1.00 

10 

. 168 

-33 

. 975 

17 

6. 16 

-108 

. 862 

-161 

.140 

-53 

.959 

15 

. 161 

-48 

.933 

22 

4.72 

-136 

. 614 

-207 

. 130 

-71 

. 890 

20 

. 155 

-64 

. 900 

26 

4. 99 

-168 

. 658 

-253 

. 132 

-85 

. 904 

25 

. 149 

-79 

. 867 

30 

4. 38 

-196 

.556 

-293 

. 127 

-97 

. 870 

30 

. 143 

-94 

. 833 

34 

3. 98 

-221 

.486 

-330 

.122 

-106 

. 835 

35 

.136 

-111 

.792 

38 

3.62 

-245 

.411 

-368 

. 114 

-123 

.781 

40 

. 128 

-127 

.742 

42 

3.23 

-273 

. 354 

-407 

.110 

-134 

.753 

-- 

— 

— 

— 

46 

2.58 

-302 

. 259 

-448 

. 100 

-146 

.685 

45 

. 117 

-144 

.679 

51 

1.92 

-328 

. 177 

-484 

.0924 

-156 

.633 

50 

. 106 

-160 

.614 

56 

1.40 

-348 

. 113 

-521 

.0808 

-173 

. 553 

60 

.0868 

-191 

. 504 

69 

. 685 

-364 

.0438 

-583 

.0640 

-219 

.438 

70 

.0688 

-220 

.400 

74 

. 820 

-380 

. 0450 

-600 

.0548 

-220 

. 375 

— 


— 

— 

78 

.917 

-400 

.0391 

-625 

.0426 

-225 

.292 

80 

.0533 

-246 

. 310 

82 

.685 

-411 

.0285 

-635 

.0416 

-224 

.285 

, 




— 

86 

.658 

-418 

.0238 

-644 

.0362 

-226 

.248 

— 


— 

— 

90 

.658 

-426 

.0201 

-695 

.0306 

-269 

.210 

90 

100 

.0414 

.0323 

-270 

-291 

. 241 
. 187 
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TABLE n. - Concluded. COMPARISON OF EXPERIMENTAL AND SIMULATION DATA FOR J85 COMBUSTOR TEST 


(e) Test E; fuel flow rate, 0. 132 kilogram per second 


Experimental data 


Frequency, 

f, 

Hz 

Spray-nozzle pres- 
sure, P„„ 

' sn 

Combustor pres- 
sure, 

Amplitude 

ratio, 

IV^snl 

Phase angle, 

Amplitude, 

N/cm 2 

Phase, 

deg 

Amplitude, 

N/cm 2 

Phase, 

deg 

deg 

7 

10. 1 

-31 

1.43 

-57 

0. 142 

-26 

12 

11.6 

-59 

1.71 

-96 

. 147 

-37 

17 

13.4 

-101 

1.92 

-146 

. 143 

-45 

22 

11.6 

-129 

1.60 

-190 

. 137 

-61 

26 

12.0 

-165 

1.66 

-239 

. 138 

-74 

30 

10.5 

-194 

1.46 

-280 

. 139 

-86 

34 

8. 89 

-219 

1.28 

-315 

. 144 

-96 

38 

7.79 

-243 

1.03 

-352 

. 133 

-109 

42 

6.78 

-266 

.807 

-388 

. 119 

-122 

46 

5.62 

-288 

.621 

-422 

. Ill 

-134 

51 

4.45 

-311 

.475 

-455 

. 107 

-144 

56 

3.50 

-331 

. 365 

-486 

. 104 

-155 

74 

2. 10 

-370 

. 168 

-595 

.0803 

-225 

78 

2. 10 

-389 

. 148 

-616 

.0707 

-227 

82 

1. 92 

-402 

. 130 

; -635 

. 0676 

-233 

86 

1. 81 

-412 

. 118 

-655 

. 0656 

-243 

90 

1.78 

-418 

. 118 

-675 

. 0667 

-256 


— 

— 

— 

1 


— 




Simulation data 


Normalized 

Frequency, 

Amplitude 

Phase angle, 

Normalized 

amplitude 

f, 

ratio, 

^V P sn' 

amplitude 

ratio, 

iV^snl/ 0 - 147 

Hz 

iV^snl 

deg 

ratio, 

IV^sJ/ 0 ' 146 

0. 966 

7 

0. 146 

-21 

1.00 

o 

o 

10 

. 144 

-30 

.986 

. 973 

15 

. 141 

-44 

.966 

. 932 

20 

. 140 

-59 

.959 

. 939 

25 

. 138 

-74 

. 945 

. 946 

30 

. 135 

-89 

.925 

.980 

35 

. 131 

-105 

. 897 

.905 
. 810 
.755 

40 

. 125 

-122 

. 856 

45 

. 116 

-138 

.794 

.728 

50 

. 106 

-154 

.726 

.708 

60 

.0865 

-185 

. 592 

. 546 

70 

.0680 

-214 

.466 

. 481 

80 

.0507 

-239 

. 347 

.460 

-- 


— 

— 

.446 

-- 


— 

— 

.454 

90 

.0375 

-261 

. 257 

— - 

100 

.0279 

-280 

. 191 


Experimental Data 


Table II contains the pertinent information obtained from the Fourier analysis of 
tests A to E. Information beyond about 100 hertz was unintelligible because of low 
signal-to-noise ratios. In addition, data around 60 hertz were obscured by high noise 
levels. 

Figures 5 to 9 contain the resultant plots of amplitude and phase for the ratio of 
combustor pressure to fuel spray-nozzle pressure. Also shown on figure 6 are the re- 
sults from the digital reduction program. All data were normalized to an amplitude 
ratio of 1. 0 at 12 hertz to eliminate any effects that rotor dynamics might have on the 
data. 

Similar response characteristics were obtained for all runs. A slower fall-off in 
amplitude was observed for the tests having higher fuel flow rates (tests D and E, 
figs. 8 and 9). The slower fall-off could be attributed to the burning zone expanding 
downstream with the increase in fuel flow. An analysis, described in the following sec- 
tion, indicates that the burning zone expanded to include about 36 percent of the combus- 
tor volume for tests D and E. For these tests, the combustor-pressure transducer 
(fig. 2) would be less affected by combustor flow and mixing dynamics. A similar phe- 
nomenon was reported in reference 5 with slower fall- offs in turbine discharge tempera- 
tures attributed to burning in the tailpipe for locally high fuel-to-air ratios. 
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Normalized amplitude ratio, Ify^srJA 132 



2 * i -300 . 

10 20 30 40 50 60 70 80 0 20 40 60 80 100 

Frequency, f, Hz 

(a) Normalized amplitude ratio. (b) Phase angle. 

Figure 5. - Response of combustor pressure to spray-nozzle pressure for test A. Fuel flow, 0. 052 kilogram per second; data normalized at 12 hertz. 





The digitally reduced data in figure 6 show apparent resonances and antiresonances 
above 60 hertz. Table II indicates that the amplitudes of both the spray-nozzle-pressure 
and combustor-pressure oscillations are down to about 10 percent of their low-frequency 
values at 60 hertz. The ratio of combustor pressure to spray-nozzle pressure was 
formed to eliminate feed system effects. The resultant amplitude ratios are subject to 
the addition of errors in measuring each pressure response. Because of these low signal 
levels, the apparent resonances and antiresonances above 60 hertz are attributed to high 
noise-to- signal ratios. Phase information seemed to be less affected by the low-signal 
levels. 



Frequency, f, Hz 

(a) Normalized amplitude ratio. (b) Phase angle. 

Figure 7. - Response of combustor pressure to spray-nozzle pressure for test C. Fuel flow, Q 092 kilogram per second; data normalized 
at 12 hertz. 
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Normalized amplitude ratio, IV^sn! /°* 146 


■30 


O 



Frequency, f, Hz 


(a) Normalized amplitude ratio. (b) Phase angle. 

Figure 8. - Response of combustor pressure to spray-nozzle pressure for test D. Fuel flow, 0. 106 kilogram per second; data normalized at 12 hertz. 



Normalized amplitude ratio, | P 4 /P sn | /0- 147 



. 4 -270 

10 20 30 40 50 60 70 80 90 100 0 20 40 60 80 100 

Frequency, f, Hz 

(a) Normalized amplitude ratio. (b) Phase angle. 

Figure 9. - Response of combustor pressure to spray-nozzle pressure for test E. Fuel flow, 0. 132 kilogram per second; data normalized at 12 hertz. 



Combustor bypass passage 



Analytical Model 

Figure 10 is a schematic representation of the system modeled on the analog com- 
puter. The annular combustor was assumed to be made up of a primary-combustion or 
burning zone and a secondary- combustion or mixing zone. Pressures and thermochemi- 
cal properties are computed at the midpoint of the burning and mixing zones (stations 4’ 
and 4, respectively). 

Fuel flow into the burning zone w^ is computed using an analog representation of 
the curve in figure 4. (All symbols are defined in appendix A. ) The fuel is assumed to 
be completely vaporized and burned in the burning zone. The flow rate of burned fuel is 
represented by w^ ^ . 

Thermochemical properties in the primary combustion zone for a range of fuel-to- 
air ratios were obtained using a digital combustion program. The resulting curves 
were implemented on the analog computer and are shown in figure 11. The portion of 
the total airflow that goes into burning w a ^ was estimated for each test using the 
manufacturers steady-state engine program and by assuming an energy balance through 
the combustor and neglecting heat transfer. By assuming that some fraction of the total 
airflow goes into burning, one can calculate the enthalpy in the burning zone from con- 
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Fuel -air ratio 


(a) Enthalpy. 


4.0xl0 3 2.4xl0 3 



(b) Gas constant. (c) Specific heat. 

Figure 11. - Thermochemical properties in burning zone as functions of fuel-to-air 
ratio. Fuel temperature, 322 K; air temperature 470 K; combustor pressure, 

30.4 newtons per square centimeter. 


ditions at the turbine entrance and compressor discharge using the standard gas mixing 
equation; that is, 


H„ 


w 


a, b 


w. 


H 3 w a. 


m 


= H„ 


w 


a, b 


w + 
a, m 


w, 


(i) 


If the calculated enthalpy H^, agrees with that obtained from figure 11 at the assumed 
mixture ratio w^/w jy then the assumed ratio is correct. This procedure yielded the 
result that the fuel-to-air ratio in the burning zone is approximately 0. 083 for all five 
tests. This value was used to estimate the steady-state airflow split. This ratio was 
later modified to match the observed zero-frequency amplitude ratios for the combustor- 
pressure to spray-nozzle-pressure responses. Analyses, described in reference 11, 
have pointed out that there is an effect of mixture ratio on the amplitude of the response 
of combustor pressure to fuel flow for rocket combustors. This result is consistent with 
the analytical modeling described in reference 12. 

Solution of the momentum and continuity equations in the air passages, burning zone, 
and mixing zone (see fig. 10) required knowledge of the size of each zone. The ratio of 


18 




burning volume V 4 , to total combustor volume (V 4 , + V 4 ) was assumed to be proportion- 
al to the square root of the steady- state ratio of burning airflow ( w a ^ + w a m j. This 
assumption was based on results obtained for representative engines using a commer- 
cial program for annular combustors (ref. 13). This assumption agrees with the gener- 
ally accepted opinion that burning takes place within about one-third of the combustor 
length. 

The burning airflow w a ^ was assumed to travel from the compressor discharge 
(station 3) to the midpoint of the burning zone (station 4’). The mixing airflow w & m 
was assumed to travel from the compressor discharge to the midpoint of the mixing zone 
(station 4). Both burning and mixing airflows were assumed to pass through a passage 
having a flow area equal to the sum of the inner and outer annuli cross-sectional areas. 
Both flow rates were computed on the analog computer using the momentum equation with 
total pressure losses proportional to the square of the flow rates; that is, 



The coefficients K, and K were adjusted for each test to match the assumed airflow 
split and experimentally measured pressures Pg and P 4 . 

The compressor discharge and combustor pressures were computed using the ideal 
gas law. At the compressor discharge, an isentropic process was assumed. This al- 
lowed, for small amplitude oscillations, the pressure to be computed from 


P 


3 ~ 



w , - w - w 
a, b a, m a, 


JdT7 + P 3 


(4) 


For the pressures in the burning and mixing zones, the time variation of the gas constant 
and temperature were included. The following equations were implemented using the 
curves in figure 11 and assuming that H = C T: 

Ir 



(5) 
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and 



( 6 ) 


Properties in the mixing zone (R 4 , H 4 , Cp^) were determined using a form of the 
mixing equation (1). A time lag for the mixing process was assumed, however, with the 
time constant calculated using average conditions in the mixing zone; that is. 


( T m S + D x 4 = £ x 3 + (1 - 0)X 4? 


v 4 p 4 . . . 

T m ” — w 4 " w a, m + w a, b + w f, b 
w„ 


(7) 

( 8 ) 


w„ 

W, 


(9) 


where X. is a gas property at station i. 

The combustor flow w CQm k was assumed to travel from the midpoint of the burn- 
ing zone to the midpoint of the mixing zone across an area equal to the combustor annu- 
lus cross-sectional area. Only momentum pressure loss was considered (i.e., P 4 = P 4 ,) 
with only a very small total pressure loss included on the analog for problem stability. 

The turbine and the bypass flow passage (see fig. 10) were assumed to be choked. 
For the turbine, the effects of the mixing-zone temperature and gas constant on flow 
were included. Including the momentum loss from the midpoint of the mixing zone to the 
turbine yielded the following expression for turbine flow: 


dr] + w^ (10) 

The effects of variations in compressor discharge temperature on the bypass flow were 
neglected and the flow rate was assumed to be proportional to the compressor discharge 
pressure. 

The compressor was modeled by a curve of airflow Wg versus pressure ratio 
P3/P2 for each value of corrected speed. Compressor speed, inlet pressure Pg, and 
inlet temperature were assumed constant for the frequency range of interest (above 7 Hz) 
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Data for the compressor curves were obtained from the steady-state engine program 
using the same compressor-face conditions (pressure, temperature, Mach number). 

Table 1(b) contains a summary of the steady- state operating conditions used in the 
analog simulation. Total airflow to the combustor was assumed to be equal to that given 
by the steady-state engine program. The listed burning and mixing airflows correspond 
to the burning-zone mixture ratios that gave the proper zero-frequency amplitude ratio 
for the combustor to spray- nozzle pressure response. Because of the high temperatures 
and potential turbine damage, no temperature measurements were made in the combus- 
tor. The listed values for enthalpy are based on the steady-state engine program and 
adjusted burning-zone mixture ratio. 

To allow flexibility in fitting the frequency response data, the following form was 
assumed for the fuel combustion dynamics: 


w. 


i, in 


-os 


f, b 2 

(a S z + bS + 1) 


(ID 


where a and b are related to the familiar natural frequency and damping ratio for a 
second- order system by 


a = 



2 ^c 


The parameters a, a, and b were adjusted on the analog for each operating condition 
to best match the experimental responses shown in figures 5 to 9. 

The mixing time constant r m could also be adjusted on the analog computer, but 
its effect was slight and good agreement was obtained using the calculated values for each 
operating condition. Appendix B contains a summary of the integral equations imple- 
mented on the analog computer. 

Sinusoidal oscillations in spray-nozzle pressure were simulated on the analog com- 
puter using a transfer -function analyzer and its associated oscillator. The analyzer was 
used to determine the amplitude ratio and relative phase angle for the response of the 
mixing-zone pressure P. to oscillations in the spray-nozzle pressure P . The 
mixing-zone pressure, rather than the burning-zone pressure P^ t , was selected because 
the experimental transducer was located in the mixing zone (see fig. 2). 
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Comparison of exper imental and computer respo nses . - Figures 12 to 16 show com- 
parisons between the experimental and simulation responses of combustor pressure 
to spray-nozzle pressure P gn - Table II contains the simulation data from which the 
plots were generated. 

In all cases, the response could be matched closely using the mixing-zone lag, as 
described in equations (7) to (9), and the form of fuel combustion dynamics given by 
equation (11). The shape of the frequency response plot was matched as well as possible 
by manipulating only the parameters a and b (see eq. (11)). The parameter a had the 
greater effect at the higher frequencies, and the parameter b was used to adjust the am- 
plitude response near the natural frequency l/2n \/a. Because of the closed- loop nature 
of the response (combustor pressure affects fuel flow through the spray nozzles), the 
dead time a had a slight effect on the amplitude ratio. However, the dead time was ad- 
justed to best match the phase angle at high frequencies for the selected values of a and 
b. Once approximate values of a, b, and a were obtained, adjustments were made to 
best fit the amplitude and phase responses. For all cases, the required fuel combustion 
natural frequency 1/27? \A/a was 55 hertz. The damping ratio b/2 yfz. had to be re- 
duced from 0. 70 to 0. 55 and 0. 50 at the higher flow rates (tests D and E, respectively). 
The slower fall-off in amplitude ratio at the higher fuel flows was attributed to the burn- 
ing zone expanding downstream (ref. 5). The need to decrease the combustion process 
damping arose from the fact that the analog simulation assumed a fixed distance (one-half 
the combustor length) between the primary combustion and pressure measurement loca- 
tions. 

As mentioned previously, the amplitude ratio at very low frequencies was matched 
by adjusting the burning-zone mixture ratio on the computer. While maintaining constant 
fuel and total airflows to the combustor, the split of airflow between the burning and mix- 
ing zones was adjusted by changing the pressure drop coefficients and K m (eqs. (2) 
to (3)). Figure 17 shows the resulting effect of the burning-zone mixture ratio on the 
zero-frequency amplitude ratio for test B. The experimentally observed zero-frequency 
amplitude ratio of 0. 158 corresponds to a burning-zone fuel-to-air ratio of 0. 0906. This 
procedure was repeated for each of the selected tests prior to the matching of the 
frequency-response data. Table in summarizes the simulation combustion parameters a, 
b, and a and the burning -zone mixture ratios that gave a satisfactory match of simula- 
tion and experimental frequency-response data. Figure 18 is a plot of the ratio of 
burning- zone volume to total combustor volume as a function of fuel flow. The volume 
ratio is based on the assumed airflow distribution (ref. 13) with the adjusted burning- 
zone mixture ratios. The dashed-curve is a fit of the data points, passing through the 
origin. Figure 18 shows the expansion of the burning-zone downstream for increased 
fuel flow to maintain mixture ratio in the range of 0. 083 to 0. 093. The higher mixture 
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Frequency, f, Hz 

(a) Normalized amplitude ratio. (b) Phase angle. 

Figure 12. - Comparison of experimental and simulation responses of combustor pressure to spray-nozzle pressure for test A. Fuel flow, 0. 052 kilogram per 
second; zero-frequency amplitude, 0.132; a fuel -to -air ratio, 0.0832; dead time, 2.60 milliseconds; mixing -zone time constant, 1.84 milliseconds; natural 
frequency, 55 hertz; damping ratio for fuel combustion dynamics, 0.70. 
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Frequency, f, Hz 

la) Normalized amplitude ratio. (b) Phase angle. 

Figure 13. - Comparison of experimental and simulation responses of combustor pressure to spray-nozzle pressure for test B. Fuel flow, 0. 76 kilogram per 
second; fuel -to -air ratio, 0,0918; dead time, 3.00 milliseconds; mixing-zone time constant, 1. 83 milliseconds; natural frequency, 55.0 hertz; clamping 
ratio for fuel combustion dynamics, 0. 70. 
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Frequency, f, Hz 

(a) Normalized amplitude ratio. (b) Phase angle. 

Figure 14. - Comparison of experimental and simulation responses of combustor pressure to spray-nozzle pressure for test C. Fuel flow, 0. 092 kilogram per 
second; zero-frequency amplitude, 0.126; fuel -to -air ratio, 0.0932; dead time, 3. 20 milliseconds; mixing-zone time constant, 1.67 milliseconds; natural 
frequency, 55 hertz; damping ratio for fuel combustion dynamics, 0.70. 
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Frequency, f, Hz 

(a) Normalized amplitude ratio. (b) Phase angle. 

Figure 15. - Comparison of experimental and simulation responses of combustor pressure to spray-nozzle pressure for test D ; fuel flow, 0, 106 kilogram per 
second; fuel -to -air ratio, 0.0863; dead time, 3. 00 milliseconds; mixing-zone time constant, 1. 37 milliseconds; natural frequency, 55 hertz; damping 
ratio for fuel combustion dynamics, 0.55. 
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(a) Normalized amplitude ratio. (b) Phase angle. 

Figure 16. - Comparison of experimental and simulation responses of combustor pressure to spray-nozzle pressure for test E. Fuel flow, 0. 132 kilogram per 
second; fuel -to -air ratio, 0.0874; dead time, 2.60 milliseconds; mixing-zone time constant, 1.33 milliseconds; natural frequency, 53 hertz; damping 
ratio for fuel combustion dynamics, 0.50. 


DO 



.40 



.101 I 1 1 I I I I 

.082 .084 . 086 . 088 . 090 . 092 . 094 . 096 

Fuel -to -air-ratio 

Figure 17. - Sensitivity of zero-frequency amplitude ratio to fuel to air ratio in 
burning zone for test B. Fuel flow, 0 l 076 kilogram per second; total combustor 
air flow, 11.3 kilograms per second (simulation data). 


TABLE m. - COMBUSTION AND MIXING DYNAMICS REQUIRED FOR MATCHING OF SIMULATION AND 
EXPERIMENTAL FREQUENCY RESPONSES OF COMBUSTION PRESSURE SPRAY NOZZLE PRESSURE 


[Natural frequency corresponding to coefficient a, 55.0 Hz. ] 


Test 

Fuel- to- air 

Steady-state amplitude 

Calculated 

Selected dead 

Selected coefficient 

Damping ratio 


ratio in 
burning 
zone 

rati< 

Experimental 

D 

Simulation 

mixing- zone 
time con- 
stant, 

T m’ 

sec 

time for best 
match of 
responses, 

sec 

for best 
respc 

a, 

2 

sec 

match of 
mses 

b, 

sec 

corresponding 
to coefficient 
b 

A 

0.0832 

0. 132 

0. 132 

1. 84xl0~ 3 

2. 60x10 3 

8. 37X10" 6 

4. 05X10' 3 

0.70 

B 

.0918 

. 158 

. 146 

1. 83 

3.00 

8. 37 

4.05 

.70 

C 

.0932 

. 126 

. 126 

1.67 

3. 20 

8. 37 

4.05 

.70 

D 

.0863 

. 146 

. 172 

1.37 

3. 00 

8. 37 

3. 18 

.55 

E 

.0874 

. 147 

. 146 

1.33 

2.60 

L _ 

8. 37 

^2.89 

.50 
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Fuel flow, kg/sec 


Figure 18. - Effect of fuel flow on burning-zone volume (based 
on match of zero -frequency amplitude ratio). 


ratios result in lower zero-frequency amplitude ratios because of the slopes of the 
thermochemical curves (fig. 11). 

Curve-fit of experimental_data. - It was demonstrated that the analog simulation 
could adequately match the experimental frequency responses. The responses obtained 
from the simulation, however, were dominated by the dead time and second-order lag 
(natural frequency of 55 Hz) used for the fuel-combustion dynamics. These values were 
assumed rather than calculated. Thus, the contribution of the simulation was to provide 
the first-order mixing dynamics and higher frequency effects due to the combustor geom- 
etry, etc. It also closed the feedback loop whereby fuel flow w^ - n was decreased as the 
combustor pressure P^ t increased. In the case of the J85 combustor operating at the 
selected test conditions, the effect of the feedback is slight. Thus, the closed-loop re- 
sponse characteristics (i. e. , natural frequency and damping) were changed only slightly 
relative to their open- loop values. Therefore, we concluded that the useful transfer 
functions could be obtained more directly by curve-fitting the experimental data. The 
zero-frequency amplitude ratio could then be obtained from the burning-zone fuel- to- air 
ratio determination. 
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Accordingly, a dead-time, first-order- lag second-order-lag combination was selec- 
ted to match the closed-loop frequency responses; that is, 


K e 
ss 


-CjS 


5 sn (C 2 S + 1) ^C 3 S 2 + C 4 S + 1 j 


(13) 


Different symbols were selected for this transfer function to distinguish it from the open- 
loop transfer function used in the simulation. 

Figure 19 shows a comparison of experimental and curve-fit responses for test B. 
Using the same procedure as outlined in the previous section, the transfer-function pa- 
rameters were ajdusted to best fit the experimental data. For simplicity, the first-order 
time constant C 2 was set equal to the calculated mixing time constant (eq. (8)): 

Table IV summarizes the results of the curve-fitting for all tests. The second-order 
natural frequency 1/2 tt \/l7c^ was 55 hertz for all tests. This result indicated that the 
combustion dynamics, as determined on the analog computer, dominate the closed-loop 
response of combustor pressure to spray-nozzle pressure. As was the case with the 
simulation damping ratio b/2 \/a, the curve-fit damping ratio C^/2 -y/C^ had to be de- 
creased at the higher fuel flows. For all tests, the magnitude of C^/2 was g rea -ter 

than b/2 \/a. Except for test E, which had the highest fuel flow rate, the curve-fit dead 
time was equal to the simulation value a. For test E, the curve-fit value was 
slightly higher (3. 0 msec compared with 2. 6 msec). 

For frequencies below about 50 hertz, the response of the J85 combustor pressure 
to fuel spray-nozzle pressure could be approximated by the dead time and a first- 
order lag having a time constant equal to the sum of C 2 and (see table IV). It should 


TABLE IV. - CURVE- FIT PARAMETERS FOR BEST FIT OF COMBUSTOR PRESSURE TO SPRAY- NOZZLE 


PRESSURE EXPERIMENT RESPONSE 


Test 

Steady-state 

Dead time, 

First-order 

Second-order lag 

Sum of 

Natural 

Damping 


amplitude 

c i, 

time 
, . a 

coefficients 

coefficients 

frequency 

ratio 


ratio 

sec 

constant, 



and C^, 

corresponding 

corresponding 




C 2’ 

C 3’ 

2 

C 4’ 

sec 

to C 3 , 

to C 4 




sec 

sec z 

sec 


Hz 


A 

0. 132 

2. 6xl0“ 3 

1. 84X10" 3 

8. 37X10* 6 

4. 97x10' 3 

6.76X10" 3 

55.0 

0. 85 

B 

. 158 

3.0 

1. 83 

8. 37 

5.21 

7.04 

55.0 

. 90 

C 

.126 

3.2 

1.67 

8. 37 

5.21 

6. 88 

55.0 

.90 

D 

. 146 

3.0 

1. 37 

8. 37 

4.74 

6. 11 

55.0 

. 82 

E 

. 147 

3.0 

1.33 

8. 37 

4.05 

5. 38 

55.0 

.70 


a Assumed equal to mixing time constant. 
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Frequency, f, Hz 

(a) Normalized amplitude ratio. (b) Phase Angle. 

Figure 19. - Comparison of experimental and curve-fit responses of combustor pressure to spray-nozzle pressue for test B. Fuel flow, 0. 076 kilogram per 
second. Curve fit values: dead time, 3.00 milliseconds,- first-order coefficient C 2 , 1.83 milliseconds; second -order coefficient C 3 , 8 . 37xl0 “ 6 second 
squared; second -order coefficient C,, 5. 21 milliseconds. Natural frequency, 55 hertz; damping ratio for curve-fit frequency response, 0.9. 



be noted that the sum of Cg and is close to the 8- millisecond lag observed for 
lower frequencies in reference 5. Below about 30 hertz, the response could be approxi- 
mated by the constant zero-frequency amplitude ratio K oc ,. 

ss 

SUMMARY OF RESULTS 

Experimental data obtained from sweep-frequency testing of a J85 turbojet engine 
were analyzed to determine the response of combustor pressure to fuel spray-nozzle 
pressure oscillations. Frequency response characteristics were displayed for all tests 
in the form of Bode plots. 

To match the experimental data and to identify the source of the observed dynamics, 
an analytical model was developed and implemented on the analog computer. The follow- 
ing results were obtained from the simulation: 

1. A dead-time - second- order- lag form for the fuel combustion dynamics was re- 
quired to match the experimental data. 

2. Dead time was relatively insensitive to fuel flow. It varied from 2. 6 to 3. 2 milli- 
seconds over the range of flows encountered. 

3. The second-order-combustion lag had a natural frequency of 55 hertz which was 
not affected by the operating conditions. The damping ratio had to be decreased for in- 
creasing fuel flow, varying from 0.70 to 0. 50. This effect was attributed to an expansion 
of the burning zone downstream toward the combustor pressure transducer. 

4. For fixed flows to the combustor, the zero-frequency amplitude ratio for the com- 
bustor pressure to spray-nozzle pressure response was a strong function of the burning 
zone fuel-to-air ratio. Burning-zone mixture ratios were computed for each test condi- 
tion. This information, together with an assumed airflow distribution into the combus- 
tor, indicated that the burning zone does expand to maintain the fuel-to-air ratio rela- 
tively constant (0. 083 to 0. 093 over the range of fuel flows considered). 

Transfer functions, obtained from curve-fitting the experimental data, were also 
determined. The following results were obtained from the curve-fitting: 

1. The response of combustor pressure to spray-nozzle pressure could be fit by a 
dead-time - first-order-lag - second-order-lag combination. 

2. The curve fit dead time matched the simulation value for all tests except the 
higher fuel flow case. 

3. The first-order-lag term was related to the mixing process in the combustor, 
and the time constant was calculated from steady-state conditions in the combustor and 
varied from 1. 33 to 1. 84 milliseconds. 

4. The second-order natural frequency was equal to the simulation fuel combustion 
natural frequency. The curve-fit damping ratio had to be decreased at the higher fuel 
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flows but exceeded the simulation combustion value for all tests. It varied from 0. 90 to 
0.70. 

5. Below about 50 hertz, the response could be fit by a dead time and a first-order 

lag. 

6. Below about 30 hertz, the response was relatively flat (+0, -2. 3 db). 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 22, 1970, 

720-03. 



APPENDIX A 


o 

A cross-sectional area, m 

a coefficient in simulation fuel- 

2 

combustion dynamics, sec 

b coefficient in simulation fuel- 

combustion dynamics, sec 

C specific heat of combustion prod- 

ucts, J/(kg)(K) 

Cj curve-fit dead time, sec 

Cg curve-fit first-order time con- 

stant, sec 

C, curve-fit second-order coeffi- 

■j 2 

cient, sec 

C 4 curve-fit second-order coeffi- 

cient, sec 

f. functional relation where 

j = 1,2, 3, 4, 5 

g gravitational acceleration, 

m/sec 2 

g c gravitational constant, 

l(kg)(m)/(N)(sec 2 ) 

H total enthalpy, J/kg 

K pressure- drop coefficient, 

(N)(sec 2 )/(m 2 )(kg 2 ) 

K choked restriction coefficient, 

(kg) (m 2 )/(N) (sec) 

K zero-frequency amplitude ratio 

s s 

K, turbine coefficient, 

(N 1 /2 )(se c)/0 t g 1 / 2 )(m 5 / 2 ) 

l length, m 


compressor speed, rpm 

2 

total pressure, N/m 
R gas constant, J/(kg)(K) 

S Laplace operator, sec”* 

T total temperature, K 

t time, sec 

3 

V volume, m 

w mass flow rate, kg/sec 

Xj fluid property at station i, 

i = 2, 3, 4, 4 T 

/3 ratio of mixing airflow to total 

combustor airflow 

y specific heat ratio 

r] variable of integration, sec 

£ c damping ratio for fuel combustion 
dynamics 

? damping ratio for curve fit of fre- 

quency response 


a 

dead time in simulation combus- 
tion dynamics, sec 

T 

time constant, sec 

w c 

natural frequency for fuel- 

combustion dynamics, rad/sec 

G 0 

S 

natural frequency for curve fit of 
frequency response, rad/sec 

Subscripts: 

a 

air 

b 

burning path from station 4' to 4 


SYMBOLS 

N 

P 


34 



c cooling path from station 3 to 4 

comb combustor path from station 4 T 
to 4 

f fuel 

in injected 

m mixing zone 


sn spray nozzle 

t turbine 

Superscripts: 

steady- state value of variable 

sinusoidally varying component 
of variable 
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APPENDIX B 


SIMULATION EQUATIONS 


Pg = Constant = P 2 


N = Constant = N 


(Bl) 

(B2) 


*a, b = t{ P 3 ~ P 4- - K b*l, bj d ” + w a, b 

\ L /b u 


(2) 



P 3 " P 4 - K m w a, m)^ + w a, m 


(3) 


P 3 = 


R 3 T 3 y 3 / * t 


v 3 '0 


f (w 2 - w a? b - ^ m - w a> c )drj + Pg 


(4) 


w 2 = f x ( N 


(B3) 


\ c = K c P 3 
in = f 2< P sn " P 4 } 


(B4) 

(B5) 


Pa.= 



C pV 4 ^ 


/ t ( w a,b +w f,b-^comb^ + P 4' 


(5) 


w 


F . = - f Fw, . (t - a) - w, . - bw, , 

f, b a Jo L f > in f ’ b f > b 


dr] 


(B6) 


*f, b = £ \ b *» * 


(B7) 
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R 4' " f 3 


w 


w 


f^b 
a, b 


w 


H 4 ,=f 4 ^ 


f,b 


w 


a, b 


c - f Wf . b 

P> 4' - f 5 

w a, b 


w 


V 4’ = 


a, b 


<*», b + "a, c>< v 4' + V 4> 


w 


^C A \ 


comb 


) A* ’ 

/comb J ° * 


P 4> d7 7 + \omb 


P 4 = 



( \omb + V m ' W t )dT ? + P 4 


R* 


tn 


ft R 3 + (1 - /3)R 4 , - R 4 drj + R 4 


H 4 = — £ f 3H 3 + (1 " $ H 4' " H 4 l dT ? + h 4 
r m u J 


-tn 


■P« r m J0 L^. + <1 -* C P«.- C Pl 


d 77 + C 


P 4 


ft = 


w 


a, m 


( *a, m + \ b + W f, b> 


(B8) 

(B9) 

(BIO) 

(B 11) 
(B12) 
(6) 

(B13) 

(B14) 

(B15) 

(B16) 
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( 10 ) 
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